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Abstract 

The atomic and electronic structures of the (001)-Si/(001)-7-Al2O3 heterointerface are inves- 
tigated by first principles total energy calculations combined with a newly developed "modified 
basin hopping" method. It is found that all interface Si atoms are 4-fold coordinated due to the 
formation of Si-0 and unexpected covalent Si-Al bonds in the new abrupt interface model. And the 
interface has perfect electronic properties in that the unpassivated interface has a large LDA band 
gap and no gap levels. These results show that it is possible to have clean semiconductor /oxide 
interfaces. 

PACS numbers: 68.35.-p,73.20.-r,71.15.Nc,02.70.Uu 
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Interfaces between semiconductors and metal oxides are playing increasingly important 
roles in advanced material science . In order to continue scaling electronic devices, a 

change from Si02 (with a dielectric constant k about 3.9) to high-/c oxides has been proposed 
for the gate dielectric in future generation metal oxide semiconductor (MOS) technologies. 
The key considerations for high-/c gate dielectrics include high dielectric constant, high 
band offsets (at least 1 eV) with respect to silicon, thermal stability, and minimization of 
electrical defects in the interface. In particular, the quality of the interface is important 
for both carrier mobility and device stability. However, control of the interface to the Si 
substrate remains a stubborn outstanding problem. For example, hafnium-based amorphous 
oxides has a bulk dielectric constant of ~ 22 4|, but its integration into the MOS gate 
stack poses substantial technological challenges 5|]. Epitaxial growth of oxides could lead 
to more abrupt oxide-Si interfaces and consequently could offer solutions for the end of the 



1) thin films have been epitaxially grown by 
Hence, AI2O3 could be a good candidate 



roadmap. Indeed, single crystal 7-AI2O3 {k ~ '. 
molecular beam epitaxy on Si(OOl) substrates 
to be used directly as a gate oxide or as a thin buffer barrier when combined with high-fc 
amorphous or epitaxial oxides. 

In another context, there have been some efforts in developing high quality crystalline 
silicon (c-Si) film on inexpensive foreign substrates such as oxides to reduce the Si material 
cost for terrestrial photovoltaic (PV) cells I?]. Previous attempts to grow single crystal Si 
on some oxides such as Ce02 failed due to the formation of Si02 |3] • Recently, Findikoglu et 
al. js] demonstrated the growth of well-oriented Si thin films with high carrier mobility on 
7-AI2O3 substrate. In addition, AI2O3 has been shown to passivate the c-Si surface efficiently 
for PV applications [Q]. These results suggest 7-AI2O3 could be a good substrate for c-Si 
solar cell growth. Therefore, detailed knowledge of the Si/Al203 interface are vital. 

Although many experimental studies have examined the growth of (001) 7-AI2O3 on the Si 
(001) surface (7-Al20a/Si), the detailed interface structure remains unclear. Theoretically, 
Boulenc and Devos [10] proposed an interface model for (001)-7-Al2O3 grown on (OOl)-Si 
surface by incorporating a defective spinel model ll| of 7-AI2O3. To obtain an interface 
without gap states, they introduced passivating O atoms to replace Si-Al and Si-Si bonds 
with Si-0 and Al-0 bonds. However, it is not clear if their proposed substrate and interface 
[lo| have the lowest total energies because only a few models were tested. It is also not clear 
if a sharp gap-states free interface can exist because the large chemical and size difference 
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between the semiconductor and the oxide. Therefore, it is desirable to obtain an improved 
microscopic understanding of the atomic and electronic structures of this important Si/Al203 
interface. 

In this Letter, we develop a new modified basin-hopping (BH) method to search for 
the lowest energy structure of the (001)-Si/(001)-7-Al2O3 interface. It is found that the 
new interface structure presents not only Si-0 bonds, but also Si-Al bonds, with all Si 
atoms 4-fold coordinated. Our density functional calculation shows that the interface is 
semiconducting with a type-I band alignment. Our results support the use of 7-AI2O3 as a 
gate oxide or a substrate for the c-Si growth. 

Our density functional theory calculations employed the frozen-core projector augmented 
wave method (PAW) [l^ encoded in the Vienna ab initio simulation package 13], and the 
local density approximation (LDA). We use a plane- wave cutoff energy of 400 eV, except 
for the search of interface structures by the BH method, where we use a soft O PAW 
pseudopotential with a cutoff energy of 212 eV. 

As a prerequisite to build the Si/7-Al203 interface, understanding the 7-Al?03 structure 
is necessary. Here, we adopt the bulk model constructed by Krokidis et al. l4] (We hereafter 



refer to this bulk model as 
defective spinel models 
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&okidis model), which has a lower energy [l5| than traditional 



IGj, and is consistent with experimental NMR and XRD results 



171]. Our test calculations also confirm the stability of the Krokidis model. The Krokidis 



model has a centrosymmetric monoclinic structure (P2i/m, No 11) with 1/4 four- and 3/4 
six-fold coordinated Al atoms. Our LDA optimization results in the following parameters: 
a = 5.479 A,b = 8.255 A, c = 7.961 A, and p = 90.645°. The lowest energy (001) 7-AI2O3 
surface [l^ based on the Krokidis model has many inequivalent surface sites (Fig. [1]). Here, 
oxygen atoms are indexed with capital letters and aluminum atoms with numbers y. It 
is noted that there is a mirror symmetry plane which is perpendicular to b and crosses Al 
3 and 4. Therefore, the O at C (D) is equivalent to the O at E (F), and the Al at 2 has 
the same environment as the Al at 5. All surface Al atoms are pentacoordinated, except 
that Al atom 1 is tetracoordinated and in a position slightly below the surface plane. All 
surface oxygens are tricoordinated if only Al neighbors within 2 A are counted. However, C 
and E oxygen atoms have an additional nearby Al atom besides the bonding Al atoms: i.e., 
the distance between Al 2 and oxygen C is 2.19 A. In this sense, C and E oxygen atoms are 
quasi-four-fold coordinated, as suggested in Fig. [T](a) by the dashed hues. Here, a four-layer 
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(not counting the tetrahedral Al atoms) symmetric slab model is adopted. After relaxation, 
oxygen D (and F), oxygen A, and oxygen B move outward from the surface by about 0.3 
A, 0.2 A, and 0.1 A, respectively. In contrast, oxygen C and E stay in the surface due to 
the strong binding with four neighboring Al atoms. It is noted that the surface is insulating 
due to the charge transfer from surface Al atoms to O dangling bonds. 

We first examine the thermodynamic stability of the interface by calculating the enthalpy 



of two possible reactions 



19|: 



|Si + AI2O3 ^ 2A1+ ISiOs, Ai7 = 2.88 eV 

1 

Si + IAI2O3 ^ |A1 + SiAlsOs, AH = 0.99 eV. 
These positive reaction enthalpies indicate that the Si/Al203 interface is thermodynamically 
stable, i.e., the formation of Si02 and silicate is unfavorable. 

The construction of the interface model is a nontrivial task. Usually, molecular dynamics 
simulations 2^ or intuition were employed for this purpose. It should be noted that molecu- 
lar dynamics simulations gives different interface structures depending on initial conditions, 
and it is almost impossible to guarantee that the constructed interface structure has the 
lowest interface energy. And it is very hard to design a good interface structure between 
two totally dissimilar materials just from chemical intuition. Therefore, we develop a new 



modified BH method 211] to determine the most stable interface structure. In conventional 
BH method, each BH run starts with a randomly chosen atomic configuration and is com- 
posed of a given number of Monte Carlo steps. In each of these, the starting configuration 
is first locally optimized to obtain an energy Ei. Then, each atom is subjected to a random 
displacement in each of its Cartesian coordinates, and a new locally optimized structure is 
obtained with energy E2. Here, Ei and E2 are the total energies from DFT calculations. 
If exp[{Ei — E2)/kBT] > r, where r is a random number between and 1 (Metropolis 
criterion), the new configuration is accepted (otherwise the old configuration is kept), and 
the process is iterated. The BH method has been widely used to search the global minimal 



structure of clusters 
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I, |2J] . However, to our best knowledge, the BH method has 
not been employed to search for the interface structure between two surfaces. 

In our newly developed modified BH method for finding lowest energy interface structures, 
we name the two slabs as "top" and "bottom", respectively (see Fig. [2]). In the case of the 
Si/Al203 interface, the Si (AI2O3) (001) slab is the top (bottom) one. The Si slab has seven 
Si layers. The top Si layer forms Si dimers and is passivated by H atoms. For the top Si slab, 



we have a rigid layer, a buffer layer, and a hopping layer. The atoms in the rigid layer can 
translate as a rigid body but the internal structure is fixed [2^. The fixed layer is held in 
place and the buffer layer is allowed to relax during the optimization. In contrast, the atoms 
of the hopping layer move as in the usual BH method but are restricted to the interface 
region. The typical value of the hopping distance of the BH simulation is about 1.5 A. Our 
test calculations indicate that swapping a Si for an Al atom is energetically unfavorable by 
about 2 eV. Thus, the bottom AI2O3 slab is divided into two parts: a fixed layer and a 
buffer layer. We note that the our modified BH method is rather general and can be used 
to search for the lowest energy structure of other interfaces. 

Considering the lateral lattice contants of the (001) Si surface [a(Si) = 5.404 A] and 
(001) 7-AI2O3 surface [a(Al203) = 5.479 A, ^(AlaOs) = 8.255 A], the best lattice matching 
is achieved by connecting the (1 x 3) Si (001) surface with the (1 x 2) (001) 7-AI2O3 surface. 
In this structure, the calculated lattice mismatch is about 1.6%. Here, the in-plane lattice 
constants of the supercell are fixed to be the theoretical lattice constants of bulk 7-AI2O3 
because 7-AI2O3 has a large Young's modulus. We perform several BH simulations for 200 
steps with different initial coordinates (the relative position between the Si surface and the 
AI2O3 surface, and the atomic positions of the atoms of the hopping layer). Finally, the 
lowest energy interface structure found from the BH simulations is refined by performing a 
full atomic relaxation of the whole system, including all atoms of the "fixed" AI2O3 layer 
and "rigid" Si layer. 

The lowest energy interface structure that we find is shown in Fig. [3l We can see that the 
dimer structure at the Si (001) surfaces is preserved as a result of the strong covalent Si-Si 
bond. We note that there is no dimer in the interface of the initial structure, while dimers 
are formed during the relaxation. At the interface, one Si atom of each dimer bonds with a 
three-fold coordinated O atom of the AI2O3 surface, whereas the other Si atom forms a bond 
with a four- fold coordinated Al atom. The Si-0 and Si-Al bond lengths are about 1.8 A and 
2.4 A, respectively. The O atoms bonded with Si move outward from the 7-AI2O3 surface 
in order to form bonds with Si atoms. We find that oxygen C and E do not bond with Si 
atoms because it is unfavorable for them to move outward due to the strong binding with the 
fourth neighboring Al atom below the surface. The binding energy between the Si surface 
and AI2O3 surface is calculated to be 2.96 eV/supercell, which indicating the strong binding 
between the two surfaces. It should be noted that there are some other nearly degenerate 
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(within 20 meV/cell) interface structures with feature similar to that shown in Fig. [31 In 
these metastable structures, other Al and O atoms are bonded with Si dimers. 

The DOS for the interface is shown in Fig. El^c). We can see that the system is semi- 
conducting with an indirect band gap of 0.46 eV. Remarkably, this value is larger than the 
LDA band gap (0.45 eV) of bulk Si. The presence of the band gap is also consistent with 
the stability of the interface. The DOS plot shows a type-I band alignment between Si 



and A 
levels 



2O3. To compute an accurate band offset, we align the energy levels using the core 



26|. The calculated valence band offset is 2.40 eV. The measured value between Si 



and a-Al203 ranges from 2.90 eV to 3.75 eV [27[. The experimental valence band offset 
between Si and 7-AI2O3 is expected to have a similar value. The discrepancy between the 
experimental result and our theoretical value are due to the different LDA error for the 
covalent Si and ionic AI2O3 but the result is qualitatively correct 28|. To gain insight into 



the electronic properties of the interface, we show the partial charge densities of the topmost 
three HOMOs and bottommost three LUMOs of the interface in Fig. [3]^a) and (b). It is 
clear that the HOMOs are mainly contributed by the directional covalent Si-Al bonds, and 
the LUMOs by the antibonding Si-0 bonds. 

It is well known that each Si atom of the symmetric Si dimer of the Si (001) surface has 
one dangling bond. On the free Si (001) surface, the tilt of the Si dimer lifts the degeneracy 
of the Si dangling bonds and a band gap opens because of the charge transfer from the 
inward Si atom to outward Si atom. In the case of the Si/Al203 interface, the band gap 
opening mechanism is totally different and much more efficient. As shown in Fig. HI the 
lone pair electrons of the surface O atom interact with the dangling bond of the nearest 
neighbor Si atom, raising the level of the dangling bond. In contrast, the high- lying empty 
Al orbital hybridizes with the dangling bond of the neighboring Si atom, lowering the energy 
level of the Si orbital. As a consequence, the Si atom bonded with the O atom transfers its 
dangling bond electron to the covalent Si-Al bond, and the interface has a large band gap. 
This bonding mechanism between Si and AI2O3 is consistent with the Bader charge analysis 
291]: the Si atom bonded with Al gains about 0.25 electrons, whereas the Si atom bonded 
with O loses about 0.40 electrons. As a result, there is some small net charge transfer (0.15 
e/Si-dimer) from Si to AI2O3. 

To investigate the kinetic stability of the interface, we calculate the energy barrier of the 
sliding of the Si surface on the AI2O3 surface. To find the transition state and energy barrier. 
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we use the "climbing image nudged elastic band" method [30|]. We consider the sliding of 
the Si surface along the b axis because the barrier of the sliding along other directions are 
expected to be larger due to the need to break all Si-0 and Si-Al bonds. The final interface 
structure is obtained from the initial structure by sliding the Si surface along the b axis by 
6(Al203)/3; the final state is almost degenerate with the initial state. In the transition state, 
there is some remaining bonding between the Si surface and AI2O3 surface: one Si-0 bond 
and two Si-Al bonds. The energy barrier of the sliding is about 2.0 eV/supercell, which 
makes the Si/Al203 interface kinetically stable. 

In conclusion, we develop a general modified BH method to search for the lowest energy 
structure of (001)-Si/7-(001)-Al2O3 interface. It is found that the interface Si dimers have 
a favorable 4-fold coordination due to the formation of not only Si-0 bonds, but also unex- 
pected covalent Si-Al bonds. Our study reveals that the Si/Al203 interface has the following 
attractive properties: (i) The interface is sharp and is semiconducting with a large LDA band 
gap; (ii) The band alignment between Si and 7-AI2O3 is type-I with both valence band offset 
and conduction band offset larger than 1.5 eV; (iii) The interface is thermodynamically and 
kineticly stable. Our results suggest that 7-AI2O3 can be used as a gate dielectric in future 
MOS technologies or a substrate for the growth of c-Si for solar cells. 

Work at NREL was supported by the U.S. Department of Energy, under Contract No. 
DE-AC36-08GO28308. 
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a 

FIG. 1: (Color online) (a) Top and (b) side view of the (1x1) (001) 7-AI2O3 surface. Oxygen (small) 
atoms are indexed with capital letters and aluminum atoms (large) are indexed with numbers. 
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FIG. 2: (Color online) The definition of various layers of the Si/Al203 interface in our modified 
BH method. 
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FIG. 3: (Color online) Interface structure and isosurface plots of the partial charge density of (a) 
the topmost three HOMOs and (b) bottommost three LUMOs of the Si/Al203 interface, (c) DOS 
plot for the Si/Al203 interface, calculated with 0.1 eV broadening. The vertical dashed line denotes 
the top of the valence band. The partial DOSs of the Si and AI2O3 surfaces are also shown. 
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FIG. 4: (Color online) Schematic illustration of the Si-Al and Si-0 bond formation and gap opening 
in the Si/Al203 interface. The valence-band maximum (VBM) and conduction-band minimum 
(CBM) of bulk Si are also shown schematically. 
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